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Abstract. The simple reflection of a light beam of finite transverse extent from a 
homogenous interface gives rise to a surprisingly large number of subtle shifts and 
deflections which can be seen as diffractive corrections to the laws of geometrical 
optics [Goos-Hanchen shifts] and manifestations of optical spin-orbit coupling [Imbert- 
Fedorov shifts], related to the spin Hall effect of light. We develop a unified linear 
algebra approach to dielectric reflection which allows for a simple calculation of all 
of these effects and lends itself to an interpretation of beam shifts as weak values 
in a classical analogue to a quantum weak measurement. We present a systematic 
study of the shifts for the whole beam and its polarization components, finding 
symmetries between input and output polarizations and predicting the existence of 
material independent shifts. 



PACS numbers: 41.20.Jb, 42.25 Gy, 42.25.Ja, 42.30.Kq 
1. Introduction 

An elementary similarity between ray and wave theories of light is the law of specular 
reflection: 'the angle of incidence equals the angle of reflection'. This rule, while valid for 
plane waves, is now understood to be only approximately true for light beams (and, in 
general, any wave beam) which has finite transverse extent. As originally discovered by 
Goos and Hanchen [1], under total internal reflection, the centre of a reflected light beam 
is slightly shifted spatially with respect to the centre of the incident beam in the plane 
of incidence. Furthermore, if the incident polarization is neither purely in the plane of 
incidence nor purely orthogonal to it, the shift has a transverse component too (Imbert- 
Federov shift) [2, 3], which may be interpreted as the spin-orbit coupling of light [4, 5]. If 
the reflection is not total, although still dielectric, the reflected beam acquires an angular 
shift, deflecting the centre slightly upon propagation, which can be both longitudinal 
and transverse [6, 7, 8]. In a closely related effect, different polarizations are reflected 
slightly differently; for incident linear polarization, different circular components are 
shifted to opposite sides of the incident plane, although the total beam centre is not 
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transversally shifted, resulting in the 'spin Hall effect of light' [9]. The separation of spin 
components of a beam is, however, not restricted to the transverse direction, and has 
also been observed within the plane of incidence [10]. Basic to all these phenomena is 
that the magnitude of the shift with respect to a suitable chosen measure is independent 
of the beam's angular spectrum, provided it is axisymmetric and tightly concentrated 
about a mean incidence angle. 

Unsurprisingly, all of these different shifts are intimately related: angular and 
spatial, in partial and total reflection, affecting the entire beam and separate polarization 
components. Our purpose here is to present a straightforward unified framework for their 
description and analysis. We derive a simple formula for the shift of any polarization 
component of the reflected beam, and show that the overall centroid of the reflected 
intensity pattern is a suitably weighted sum of the shifts of orthogonal polarizations of 
the final beam. This form is naturally expressed as a complex vector, and writing it 
in this way reveals a close relationship with the quantum-mechanical notion of weak 
values and weak measurements of operators [11, 12]. In this approach, the effect of an 
operator on an initial state is not considered overlapped with the same initial state, as 
with a regular quantum expectation value, but rather with a particular, postselected 
state which is different from but not orthogonal to the initial state (and neither the 
initial nor the final state is an eigenstate of the operator). This has been previously 
discussed and experimentally measured for the 'spin Hall effect of light' [9], where a 
weak measurement was used to show that beam components undergo shifts that the 
overall intensity centroid does not experience. 

In our case, the initial state preparation corresponds to a uniform initial polarization 
E, a postselecting measurement state from a polarization analyzer F (in general 
unrelated to E), and the operator in question is the reflection matrix, which acts on 
different plane wave components of the beam differently, since the angle and plane 
of incidence depend on the direction of the plane wave component. For a narrow 
paraxial beam, this dependence is so small that the polarization couples only weakly 
to space, so the spatial beam shift is analogous to shift of the pointer. The weak value 
of an operator is in general complex, and it has been proved that the imaginary part 
corresponds to a shift in the conjugate variable [13, 14]. If the pointer corresponds to 
the real space representation of the beam, then its conjugate variable is the Fourier 
space: the imaginary part of the weak value is the angular shift. The full extent of the 
connection between weak values and beam shifts extents even further and deserves its 
own exposition in a more general framework, which we present in a companion paper 
[15]. 

Here, we consider explicitly the resulting shifts with input and analyzer polarizers 
(E and F) chosen from three natural bases: the transverse electric (TE) and magnetic 
(TM) ('+'), circular ('O') an d 45° and 135° linear ('x') bases, and we find a symmetry 
in the shifts between input and output polarizations common to total and partial 
reflection. Surprisingly for a phenomenon which can be used for optical biosensing and 
nano-probing [16, 17] we also find shifts of polarization components that are material 
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independent. Furthermore, our unified formula shows that, when measured in the 
appropriate units, the formula for the angular shift is independent of the shape of 
the spectrum of the incident beam. Our approach is self-contained, and our derivation 
clarifies aspects of previous descriptions of beam shifts. 

The physics of beam shifts depends wholly on the Fresnel reflection coefficients r, 
which are the scalar weightings an electromagnetic plane wave acquires on reflection 
from a planar reflection surface. We assume throughout that the reflection is dielectric 
(i.e. unit relative magnetic permeability ji and real relative permittivity e), and an angle 
of incidence 9 > 0. In general, there are two eigenpolarizations for plane wave reflection, 
's' corresponding to (electric) polarization perpendicular to the plane of incidence, and 
'p' for (electric) polarization in the incidence plane, with two different Fresnel coefficients 
[18]: 

cos 9 — i-\/ sin 2 9 — n 2 cos 9 — y ' n 2 — sin 2 9 



cos 9 + i a/ sin 2 9 — n 2 cos 9 + \/ n 2 — sin 2 9 



(1.1) 



n 2 cos9 — ivsin 2 9 — n 2 n 2 cos 9 — y n 2 — sin 2 9 
r p = = (1.2) 

n 2 cos 9 + i v sin 2 9 — n 2 n 2 cos 9 + y n 2 — sin 2 9 
Here, n is the relative refractive index (the refractive index of the medium in which 
the light propagates is in the denominator); that is n > 1 for reflection by a denser 
medium. The reflection coefficients are each given in two forms, each useful depending 
on whether sin 2 9 — n 2 ^ 0. The first assumes that sin 2 9 > n 2 , appropriate for total 
internal reflection (where 9 is greater than the critical angle arcsinn), in which case 
the reflection coefficient is complex and unimodular, |r| = 1. In the other case, partial 
reflection (which also has some transmission, which not considered here), sin# < n, so r 
is real. Although the formulaeare similar, r p = at the Brewster angle arctan n, whereas 
r s 7^ for dielectric reflection. 

The spatial and angular Goos-Hanchen (GH) shift within the plane of incidence can 
be intuitively understood from a two-dimensional picture in which s and p polarization 
can be treated separately. Any two-dimensional wave packet or beam consists of a 
spectrum of plane waves with different incidence angles, so each acquires a different 
reflection coefficient (for each s and p component of the polarization). When the beam 
is tightly concentrated around a common incidence angle 9q, ai l beam shifts are a first- 
order effect in the spectral width of the beam, and their nature (spatial or angular) 
depends on how the polarization-weighted sum of r s and r p varies over the spectrum: a 
variation of the imaginary part alone leads to a spatial shift, and real variation, to an 
angular shift (see below, Fig. 1). Because both spatial and angular Goos-Hanchen shifts 
are usually small when compared to a generic measure, they can be seen as corrections 
to the laws of geometrical optics [7, 19]. 

To understand the Imbert-Federov (IF) shift, which is transverse to the incidence 
plane, it is necessary to use a vectorial, three-dimensional picture, and it is crucial to 
distinguish between the global direction of the beam's electrical field and the polarization 
of the individual plane waves in the angular spectrum. Owing to transversality, the 




Figure 1. Schematic illustration of 'virtual reflection' and spatial and angular 
Goos-Hanchen shifts, a) The angular spectrum of the incident light beam (green) 
is schematically shown as colour gradient. The 'virtual reflected beam' (red) is formed 
by the geometric reflection with respect to the interface (x-axis) of each plane wave 
in the spectrum, b) For a unimodular reflection coefficient (1.1), (1-2) the physically 
reflected beam appears shifted by distance D parallel to the propagation axis of the 
virtual reflected beam, c) For a real reflection coefficient the physically reflected beam 
is deflected by an angle A due to the modulation of the amplitude for each plane wave 
in the spectrum upon reflection. 



electrical field for each plane wave in the spectrum has a different direction, and 
therefore different s and p polarization components with respect to their local plane 
of incidence. This is the origin for the 'spin-orbit coupling of light' [5] and is responsible 
for spatial and angular Imbert-Fedorov shifts [2] . In contrast to the in-plane shifts, these 
transverse shifts do not depend on the change of the reflection coefficient with the angle 
of incidence, but rather on the difference between the reflection coefficients for the s 
and p components. 

In our calculation, both effects appear in a unified framework by virtue of a general 
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scattering or reflection dyadic R, which projects the electric field of each plane wave 
component into its local s and p direction, and applies the appropriate Fresnel coefficient 
to the reflected wave. The main formulae for the total beam and component shifts are 
given together in Section 2, and are interpreted in terms of weak values of the beam. 
These expressions are derived in Section 3. 

We compare the total and component shifts for three naturally incompatible bases 
of polarizations for the incident and reflected beam in Section 4: the + basis given 
by the transverse electric and magnetic polarizations (i.e. the eigenpolarizations of a 
plane wave on the beam axis), the O basis of circular polarizations, and the x basis 
of linear polarizations diagonal to the + basis. This choice of three bases (represented 
by three orthogonal directions on the Poincare sphere [18]) is clearly special, since 
each polarization is made up of equal weightings of each vector in the other two bases: 
for each one of these polarization states, a polarization measurement in a different 
basis cannot discriminate between this vector and its orthogonal partner. There are 
additional symmetries that occur between different polarizations, some of which are 
specific to whether the reflection is total or partial. 

2. Formulae for total and component beam shifts, and interpretation as 
values 

For ease of reference, we present the main formulae for all of the various shifts, including 
the well-known total beam shift, and the newly derived scalar component shift. A 
self-contained calculation for all of them is given in Section 3. 

We think of two different setups, which lead to two different types of shift: the 
first is the 'centre of mass' shift of the entire vector beam, that is without specifying 
a polarization analyser for the reflected beam. The second is the shift of a polarized 
scalar component of the beam, as filtered out from the non-uniformly polarized reflected 
beam by use of an analyser. The fact that different polarization components undergo 
different shifts has been called the 'spin Hall effect of light' [20] when these shifts are 
transverse, although as we will describe for the general case, both longitudinal and 
transverse relative shifts between different polarization components are typical. The 
total centre of mass shift is, in fact, an average between every pair of orthogonal 
polarization components, weighted by the intensity of that component (as we will show 
explicitly later). 

We use the coordinates of the reflected beam (as below, Fig. 2), with x-, y- and z- 
directions corresponding to the axial TM (p-direction of the axis), TE (axial s-direction), 
and propagation directions respectively. (We use TE and TM to avoid confusion with s 
and p, which are defined differently for each plane wave component.) The polarization 
of the incident beam in the TM-TE plane is written E and an output polarization 
component is determined by the analyzer polarization F which in general are complex 
for elliptic polarization, and given in terms of TE, TM components, 

E = (E m ,E e ), F — (F m , F e ). (2.1) 
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Our formulae will be written as a complex dimensionless vector T>, whose 
corresponding spatial shift D and angular shift A are given by 

D = -hmT>, A = (5 2 )ReT>, (2.2) 
k 

where (<5 2 ) is the second moment of the spectrum in Fourier space, and k is the optical 
wavenumber. The ^-component of these two real vectors, in the TM direction of the 
beam, is therefore the longitudinal GH shift, and the second, y-, component is the 
transverse IF shift. The shift formulaedepend on the beam axis incidence angle 6q, the 
mean reflection coefficients r p and f s evaluated at this angle, and their derivatives f' p , 
f' s , with respect to 8 at 6q. The overall beam shift depends on the incident polarization 
Eq, and the shift of a single polarization component depends on E and F . 

The expressions for the shifts can be written more compactly with the definition of 
vectorial quantities relevant to reflection. The 'mean reflection matrix' in the TM-TE 
basis is 

*-("?".) (2 ' 3) 

and the first-order term in the small incidence angle expansion of the full matrix depends 
on two matrices, 

R<ui=( ' p I, ), RiF = -cottf | ,_ ° {fs+ n fp) I- (2-D 




(f s + f p ) 

The diagonal R G h is responsible for the longitudinal shift, whereas the off-diagonal R IF 
determines the transverse shift. 

The shift to the total beam T> t is simply 

t> _ ( ■ R^Rgh -EE*- R*R IF ■ e\ 

V E R R E E* ■ R R E J 

where R^ denotes the conjugate transpose of the mean reflection matrix R and E* is 
the complex conjugate of the polarization vector E. In the derivation of this expression 
in Section 3 we write T> t slightly differently, but the expression above is equivalent to 
(3.20) using the convention of the 2-component shift vector in which the first and second 
components refer respectively to the shifts in the longitudinal and transverse direction. 

The shift for a polarization component T> c is given by a similar expression, but now 
depending on the output polarization F, 

/ F*.Ran-E F*.R 1F .E \ 

V F* R E F* R E J y ' 

Both shift expressions could be further simplified by writing a 2- vector of 2 x 2 matrices 
(Rgh, Rif), but we will keep the separation explicit in the following. 

The form of equation (2.6) is very similar to that of a the 'weak value' of a quantum- 
mechanical operator A with initial state \E) and final state |F). The usual definition of 
the weak value of A is [11]: 
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The significance of this quantity is that it represents the outcome of the operation A\E), 
postselected into a new state \F) rather than compared with the initial state \E). The 
weak value is of course only defined when \F) is not orthogonal to \E) (otherwise the 
denominator would blow up the whole expression), but can take on arbitrarily large 
values when \F) and \E) are almost orthogonal, as the overlap in the numerator is not 
necessarily particularly small, although the denominator can approach zero arbitrarily. 
Hence weak values can be arbitrarily large even when the spectrum of A is bounded, 
and as such are sometimes called 'superweak' [21]. It is crucial that the interaction of 
the initial state with the rest of the system is weak, so the final system A \E) is largely 
undisturbed by the action of A; this weak interaction should involve changing the state 
of a different, pointer state whose motion indicates the weak value, given by KeA w . We 
will discuss this connection briefly below; more details of the connection are given in 
the companion paper to this one, [15]. 

We consider a well collimated beam with a narrow axisymmetric spectrum a (5) 
which is parametrized by a single angle 5, which denotes the angle of every plane wave 
in the spectrum with the central wave vector or the beam axis (see Fig. 2). This beam 
is in an initial state, namely a wavepacket with a homogeneous initial polarization E, 
represented by a single point in polarization space (the Poincare sphere [18]). After 
reflection, the beam is strongly localized in polarization space, but is no longer a single 
point (and, when E is not an eigenpolarization, the state out is not the state in). 
Since the initial beam was concentrated in a particular initial direction, the spread of 
polarization in the reflected beam is small, but postselecting by a particular component 
F can reveal a larger effect than is present in the overall beam (the total shift). 

Furthermore, since the spread in Fourier space is narrow, the beam width is 
large, much larger than the small (wavelength-scale) shift in the beam's centre. The 
configuration space position of the beam thus corresponds to the measurement system, 
interacting weakly (due to the narrow beam) on reflection with the polarization degrees 
of freedom: ReiX> c corresponds to KeA w . 

The imaginary part of the weak value (2.7) is related to a shift in the canonically 
conjugate variable corresponding to the pointer, i.e. the direction of the beam (position 
in Fourier space). In Ref. [14], it is proved that, under a quantum weak interaction, 
the weak shift in direction should be proportional to Im(A w )Varfc, where Var^ is the 
variance of the wavenumber k, which — in our case of an axisymmetric spectrum - 
is given by (<5 2 ), the second moment of the spectrum in the Fourier space. The case of 
the coupling of direction and polarization in beam reflection is similar to von Neumann 
measurement: interaction is instantaneous in propagation distance (i.e. time), and the 
weak interaction Hamiltonian (the reflection operator) depends linearly on direction k. 
Indeed, not only is the angular shift proportional to the imaginary part of iX> c , but the 
proportionality factor is precisely the variance of the spectral distribution in Fourier 
space. 

In their experimental verification of the spin Hall effect of light, Hosten and Kwiat 
[9] interpreted the very large spatial shift on transmission of an almost vertical analyzer 
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state with a horizontal initial state, on transmission through a dielectric object (we 
make the reasonable assumption that there are similar transmission counterparts to our 
equations in Section 2). In particular, the sign of the shift was seen to depend on the 
sign of the small perturbation angle from vertical, and of course the shift was not defined 
when the analyzer was truly vertical. 

Of course, the denominators of (2.6) and (2.7) differ, which in the beam shift case 
depends on the mean reflection matrix R, which is (possibly complex) diagonal in the 
TM-TE basis. This corresponds to the fact that, unless the initial beam is polarized in 
the TE or TM direction, the disturbance to the polarization on reflection is not small, 
and is given by R. The shift is therefore weak with respect to the system (without 
interaction) evolving into \E) = H\E), corresponding to the polarization of the on- axis 
plane wave polarization after reflection R • E. With this replacement, we can see that 
the beam shift truly corresponds to the weak value of the operator we call the 'Artmann 
operator' 

IAH->) (F\AH-'\E) (F\6\E) 

preselected by the result of the noninteracting evolution operator H acting on \E), and 
postselected by |F). Note that the notation here differs slightly from [15]. 

Superweak values are therefore guaranteed when the postselection analyzer F is 
almost orthogonal to R ■ E, rather than almost orthogonal to E itself. For the choice 
in [9] , E is an eigenstate of R, so in this case only should F be almost perpendicular 
to E. 

A different, arbitrarily large shift is possible, distinct from the aforementioned 
'superweak shift': the mean reflection matrix R itself may be close to singular, which 
occurs when 9q is close to the Brewster angle (at which f p — 0). To take advantage of 
this large shift, either the analyzer or the numerator should be wholly longitudinal, say 
the analyzer F. In this case, F* ■ R ■ E = —f p E m , which is arbitrarily small regardless 
of the initial polarization. Of course, this discussion is only valid in the first order 
in 5 approximation we have developed, and more sophisticated mathematics has been 
employed to describe reflection at Brewster incidence [22, 23, 24]. 

Therefore, given the presence of the mean reflection operator in the denominator 
of (2.6), there should be a finite but nonzero shift between any pair of orthogonal 
polarizations. After the derivation of the shift formulae, given in the next section, we 
will discuss the spatial and angular shifts that occur between states in the +, O an d 
x bases, for both total and partial reflection. However, first we give a derivation of the 
beam shift formulae. Readers uninterested in these derivations, but are interested in 
their consequences, may skip to Section 4. 



3. Calculation of beam shift formulae 



Our approach to the calculation of the shift formulaegiven in the previous section is based 
on the introduction of a 'virtual reflected beam', which is in reflected beam coordinates, 
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constructed by 'mirroring' every plane in the angular spectrum as in geometrical optics 
without taking the reflection coefficient into account. Effectively, the virtual beam is 
the reflection of the incident beam with r s — l,r p — 1 for all directions, but otherwise 
the virtual beam incorporates the direction change of each plane wave component on 
reflection. It is therefore the virtual beam which is depicted in Fig. 1 (a). The real 
reflected beam follows from applying the reflection coefficients (1.1), (1.2) to each plane 
wave component of the virtual beam separately. 

We find the virtual beam a useful conceptual object in considering beam shifts, 
particularly the angular shift. The origin of the angular shift comes from the fact 
that different plane wave components acquire different amplitude weightings, as the 
real reflection coefficient depends on angle. In Fourier space, the shift comes from a 
real gradient applied to the spectrum, as in Fig. 1(c), rather than a translation of the 
overall direction (which would be closer to the spatial shift). The virtual beam is the 
same, independent of the physics of the different physical regimes considered (spatial 
or angular, total or partial reflection). We note that the virtual beam utilized here 
is slightly different from the one we describe in [15] - there, the virtual beam has a 
different overall polarization from that considered here. 

We therefore work in beam coordinates defined above, with x, y and z the TM, TE 
and propagation directions. In coordinates based on the dielectric interface (implicit 
in Fig. 1), the reflection angle of the beam is 9q, which therefore corresponds to the 
^-direction. In spherical polar coordinates based around the propagation direction, we 
define an azimuth angle a and colatitude 5, which is zero on the beam axis. Each plane 
wave making up the beam has a wavevector of magnitude k, and is parametrized by a 
and 5, 



The spectrum, giving the real amplitude weighting of the plane waves, is denoted cr(5), 
which is concentrated around 5 = 0, that is, the vector (0,0, 1) along the propagation 
axis. We also assume that the spectrum is independent of a (i.e. axisymmetric), and 
is square-integrable, but apart from this (and the concentration around 5 = 0), the 
spectrum is otherwise completely general. By transversality, the initial polarization of 
the axial (0,0, k) component beam must be perpendicular, with components E m ,E s in 
the TM, TE directions (and given by the 2- vector E in this plane). Since the incident 
beam has a constant overall polarization given by the constants E m , E e , the polarization 
of each plane wave component must have these components in the TM, TE direction. 

The polarization of the incident beam, and hence the virtual beam, is uniform 
polarized. However, since the plane of incidence varies around the beam azimuth a, 
the precise resolution into each plane wave's particular s and p direction varies around 
the beam (Figure 2). A direct consequence of this variation is that the reflected beam 
will no longer be uniformly polarized, which is why choosing different settings for the 
polarization analyzer F can lead to very different beam shifts. The shifts themselves 
only depend on terms to first order in 5 in the reflection coefficient and local incidence 




(3.1) 




Figure 2. Schematic explanation of the beam coordinate system x, y and z, where 
z points along the propagation direction of the beam, x lies in the global plane of 
incidence and y is transverse to both x and z. a) shows the cone geometry used to 
parametrize the angular spectrum by 8 and a and illustrates that the local angle of 
incidence 9 varies around the angular spectrum, b) The local planes of incidence are 
not coplanar with the x, z plane and this is why the decomposition of the polarisation 
into the local s and p directions varies around the spectrum. 

direction; higher-order terms have a smaller effect which we ignore. 

One such ignorable effect comes from the fact that each plane wave acquires 
a component proportional to tan 5 in the (0,0,1) direction, which couples into the 
transverse plane by projection into the local s,p directions: this coupling is also of order 
5, so this extra spin-orbit effect plays no significant role in the transverse shifts, though 
it affects the axial polarization [25]. The whole analysis can therefore be restricted to 
the two transverse dimensions alone. The virtual reflected beam in real space, E v , is 
therefore given by 

E v (r) — / dS / da a(5) sin 5 exp[i/c sin 5(x cos a + y sin a)]E, (3.2) 

J0 J -IT 

where x,y are TM, TE beam coordinates, and we ignore z-dependence. We are not 
concerned with the upper limit in the S integral, as it is cut off by the narrow spectrum 
a(5). 

In beam coordinates, the reflection plane is now oblique, as depicted for a cone of 
wavevectors with constant 5 in Fig. 2. The intersection of the /c-space cone and the 
plane is thus a conic section, and the precise angle of incidence of each vector on the 
cone varies with a. If the reflection angle of the beam axis is 8q, the normal vector to 



Generalized shifts for polarization components of reflected light beams 



11 



the reflection plane in beam coordinates is 

n = (— sin# , 0, cos# ) (3-3) 

so the beam coordinate ^/-direction is perpendicular to the plane of incidence. The angle 
of incidence of each vector on the cone is therefore given by cos 9 = n ■ k or 

cos #(#o ; a ) — cos $ cos $o — cos a sm $ sm $o- (3-4) 

Whatever the overall incoming polarization E, it is resolved into the local s,p 
directions for each component wavevector labelled by a, 5, each of which also acquires the 
appropriate reflection coefficient r s ,r p from (1.1), (1.2), with its local angle of incidence 
provided by (3.4). As we are only interested in effects linear in 5, these coefficients 
effectively become, for j = s,p, 



r 



j ~ Tj + 5cosarj ss r,- exp(5 cos a r'j/rj). (3.5) 

The second approximation is only valid if fj ^ 0, i.e. not at the Brewster angle. 

For a given wavevector k(S, a), the s polarization perpendicular to the local plane 
of incidence is given by the direction of the vector product with the normal direction n, 
namely 

nxK(5,a) 

e s (5,a) = - — — . 3.6 

\n x k(o, a)\ 

The p polarization vector is the cross product of e s with k, 

e p (5, a) = e s (5, a) x n(5, a). (3.7) 

These vectors give rise to projection matrices in the TM,TE plane, which are then also 
approximated up to order 5, 

, s ( —5 cot 6 sin a \ , n oX 

P s = (e s ®e s w . , 3.8 

\ — o cot0 sma 1 / 

, x / 1 5 cot #o sin a \ , n n . 

P p = {e p ®e p w I . 3.9 

\ o cot 6*o sina / 

Clearly, the variation of order 5 in these projection matrices is proportional to sin a, 
as can be seen in Fig. 2. When a = 0, or ir, the local plane of incidence agrees with 
the plane of incidence of the beam axis, so there is no additional polarization effect. 
This dependence on sin a has been interpreted as a spin-orbit effect [5], as it implicitly 
involves considering the local polarization of the beam (spin) around the beam azimuth 
(orbit). 

The total transverse reflection dyadic, then, is the sum of the projection matrices 
weighted by the appropriate reflection coefficient. Before we do so, however, we point 
out a subtlety made explicit by the introduction of the virtual beam: since the virtual 
reflected beam comes from reflection, it involves a change in the handedness of the 
coordinate system (right-handed circular light is reflected as left-handed, and vice versa). 
Nevertheless, it is natural to make the coordinate systems attached to the incident, 
virtual and real reflected beams right-handed, so for every plane wave, the triple of 
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wavevector fe, electrical field, and magnetic field has a positive scalar triple product. 
To account for this discrepancy we choose introduce an additional sign change in the 
reflection coefficient r p , which leaves the transverse TE (y) direction unchanged. With 
this sign convention the total reflection dyadic is 

/ -f p \ / -f'cosa -(r s + r p ) cot 6» sin a 

R = r s P s -r p P p ^ " _ \+S\ ( _ \ 

\ r s I \ — (r s + r p ) cot # sina r s cosa 

= R + 5cosaR G H + 8 sinaRip, (3.10) 

plus terms of order 5 2 and above. The second line provides the origin of the matrices 
introduced in (2.3), (2.4) above. 

We are now in a position to compute the beam shifts given in Section 2 above. The 
total reflected beam E r (r) is constructed by putting the first-order polarization dyadic 
(3.10) into the virtual beam (3.2), 

E r (r) — / d5 / da sin5o-(5) exp [i/c sin 5 (a; cos a + y sina)]R • E. (3-11) 

JO J-TT 

A scalar component of the reflected field is selected by projecting this beam onto the 
direction of a polarization analyzer F (whose components are constant). The reflected 
scalar component is thus given by 

F*-E r (r)= [ d5 [ da sin 5a(5) expfiA; sin 5(x cosa + y sin a)]F* RE. (3.12) 

JO J-TT 

The scalar F* ■ R • E thus determines the nature of the general shift and it is not 
necessary to calculate the integral in Eq. (3.12) explicitly. 

The scalar component of the reflected beam (3.12) only depends on an angle- 
dependent reflection scalar 

R = F* ■ R- E w R + 5(cosaF* ■ R GH • E + sinaF* ■ R IF • E), 

= R[l + 5T> C ■ (cos a, sin a)] , (3.13) 

where D c is the complex shift vector dependent on F, E, defined in (2.6), and 
R = F* ■ R • E, the mean reflection coefficient for the scalar component given by F. 
Using the same reasoning as in the expansion of the reflection coefficient in Eq. (3.5), 
R may be expressed as an exponential, 

R w (F* ■ R • E) exp [5T> C ■ (cos a, sin a)] . (3. 14) 

The trigonometric factors cos a and sina therefore determine which terms in the 
exponential are responsible for the respective shifts in the transverse and longitudinal 
directions, because inside the exponential in Eq. (3.12) cos a couples to x whereas sina 
couples to y. Alternatively, the ^-component of T> c contains all those terms in the 
exponent of R in Eq. (3.14) which have cos a as a common factor, while the y component 
collects all the terms which multiply sin a. 

To see that this complex vector really gives rise to both angular and spatial shifts 
it is easiest to split T> c into real and imaginary parts explicitly, 

R ~ R [1 + 5Re(T> c ) ■ (cos a, sin a)] exp [i5Im(X> c ) • (cos a, sin a)] , (3.15) 
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with the imaginary part in the exponent. Substituting this expression for F* HE into 
the expression (3.12) of the reflected beam and identifying 5 ~ sin<5 in the exponential, 
we construe the square brackets containing the real part of T> c as a direction-dependent 
modification of the spectrum a (5), and the exponential with the imaginary part as a 
shift of the spatial variables within the Fourier kernel. Setting D c = (D cx ,D cy ) (see 
Eq. (2.2)) in the transverse plane gives 

F* ■ E r (r) — R fdS [ da sm5a(5) [1 + <5Re(X> c ).(cos a, sin a)] 

Jo J-n 

x exp[i/csin<5 (r — D c ) ■ (cosa, sin a)]. (3.16) 

The effect of the imaginary part of T> c has therefore been to translate the transverse 
spatial coordinates of the beam by D c , given in (2.2) for the particular choice of 
component. The change in the spectrum does not affect the spatial variables to leading 
order in 8. 

The angular shift is calculated as the expectation value of the spectrum of the 
reflected component as a distribution in Fourier space. This distribution is the modulus 
squared of the spectrum of (3.16), which is equal to \o~(8) | 2 , unaffected by the unimodular 
spatial shift, which only affects the spectrum's phase. We therefore calculate the centroid 
of an angle vector (S x , 5 y ) = 5(cos a, sin a), whose components represent the longitudinal 
and transverse coordinates in Fourier space about the mean propagation direction (0, 0), 
which is also the centroid of the virtual reflected beam. The correctly-normalized 
expectation value in Fourier space is therefore 

J d5 f^da |cr(<5) | 2 sin 5 [1 + 5Re(T> c ) ■ (cos a, sina)] 2 5(coso;, sin a) ^ 

f Q d5 da \cr(S) | 2 sin 5 [1 + 5Re(T> c ) ■ (cos a, sin a)] 2 

For the evaluation of this integral the term [1 + 5Re(T> c ) • (cos a, sin a)] 2 is crucial, and 
keeping only the the two lowest orders in 5, it can be approximated to 1 + 2<5Re(X> c ) ■ 
(cos a, sin a) both in the numerator and denominator. Utilising the fact that cos and 
sin are even and odd functions it is straight forward to evaluate the integral which gives 
the final centroid as 

(6 X ,6 V )=A C , (3.18) 

defined for the beam's component by (2.2), with the second spectral moment (5 2 ) given 
by 

/r9v Ld5\a(5)\ 2 5 2 sin 5 , x 

As indicated above, using the second spectral moment as the appropriate unit for the 
angular shift makes the the angular shift independent of the spectrum when the shift 
is measured with respect to (S 2 ). This measure is different for a two dimensional wave 
packet, accounting for the factor 2 in the formula for angular shift of a quasi-scalar 
'sheet' of light [26]. 

With the shifts now defined for any polarization component of the output beam, 
it is possible to rederive the shift of the overall intensity pattern of the vector reflected 
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beam E r . The complex shift vector X> for this total 'centre of mass' shift is a sum of the 
shift vectors for two orthogonal polarization components as analyser settings, weighted 
by intensity. If one such analyzer polarization is given by F = (F m ,F e ), its orthogonal 
partner can be denoted F± = (F*, —F^), so F* -F± = 0. The mean reflection coefficient 
R = F* H E is, to leading order, the amplitude of the scalar component of the reflected 
beam corresponding to F, and 

R ± = F* ± • R • E = -F m E e f s - F e E m f p (3.20) 

the amplitude for the orthogonal analyzer component F±, which has a complex shift 
vector T> c _|_ given by (2.6) with Fj_ instead of F. The sum of these two complex shift 
vectors, weighted by their respective intensities |-R| 2 , \R ± \ 2 is therefore the total reflected 
beam is given by 



T> = 



\R\ 


2- 


T> c +\ 


R±\ 






R\ 


2 + 


\R± 


2 



(3.21) 



which gives the F, .F^-independent form (2.5) after straightforward but tedious 
algebraic manipulation. We have therefore given an alternative definition of the total 
shift vector of a vector beam as a weighted sum of its more fundamental scalar 
components. 



4. Centre of mass and scalar component shifts in the '+ O x ' bases 

In this section we illustrate the variety of beam shifts present in the equations (2.6) 
and (2.5) for different choices of incident and analyser polarization. We choose three 
bases of orthogonal polarizations as settings for the incident polarisation E and the 
analyser polarisation F and consider all possible combinations of pre- and postselection 
for these bases. The three bases are TM and TE polarisation ('+'), right and left circular 
polarization ('O') an d the two diagonal polarizations of 45° and 135° degrees ('x'). On 
using the same notation for the components of the polarisation above, the '+0 x ' bases 
are explicitly given by the following polarization vectors: 



E*+ = F*+=l , E° = F° = , eS = fS = -L[ (4.1) 






E t = F t = [ - ] ^ E = F° = ±y* i j,Er* = F** = ^=[ \ ) (4.2) 

An algebraic representation of the different beam shifts as 2D vectors for every 
combination of polarizer and analyser setting in the '+ O x ' bases is given in Table 1. 
The table parts naturally in four blocks: the eigenpolarization block ('+/+'), in which 
both polarizer and analyser are set to either TE or TM, two blocks in which either 
the polariser or the analyser is set to an eigenpolarization ('+/ O x ') or ('+/ O x ')> 
and one block which contains the combinations of circular and diagonal polarizations 
CO x / O x ')- Because of the change of handedness on reflection we have switched left 
and right circular polarisation as analyser settings in the last block, which makes the 
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Table 1. Table showing in an algebraic form the resulting complex shift vectors 
for the total shift in Eq. (2.5) and and shift of a scalar component Eq. (2.6) 
for different input polarizations E and in the case of the component shifts also 
different analyser polarizations F. The formalae for the expressions d p , d s as well 
as p, q 7 r, u, v, a, c, e, g, f, h are given in Table 2. 



Analyser 
setting F 



TM 



Polarizer setting E 


























TE 



RC 



LC 



45° 



135° 



no analyser 



d s 




I 



TM 



TE 






d s 




iu 

d s 
iv 



—iu 

d s 
—iv 



p 
u 

d s 

—v 



up 

— U 

d s 
v 



LC 

C RC 
45° 



I 



135° 



ill 



-iu 



d, 



"p 
-u 



d s 

IV 

ds 
—iv 

ds 

—v 

d s 

V 



-ib 

9 
h 



e 
f 

9 
h 

c 

-d s 

a 




a 


c 
ds 



whole table symmetric. Symmetries along the rows of the table are therefore reflected 
in symmetries along the columns and it is often convenient to discuss only one. 

The table follows the notation in this paper in that the first component of each 
vector refers to the shift in the plane of incidence (longitudinal) and the second 
component to the shift out of the plane of incidence (transverse). As before the 
real part of each component determines the angular shift, while the imaginary part 
determines the spatial shift and this is why we may expect to find symmetries in the 
spatial shifts involving diagonal polarisation also mirrored in the angular shifts for the 
circular polarization. For comparison the shift of the total beam, without filtering out 
a scalar component by the use of an analyser, but for different incident polarizations, is 
shown in the top row. 

The entries in Table 1 are derived solely from the linear algebra in section 2 and 
are valid for both total and partial reflection. The most simple block to explain is the 
eigenpolarization block where the only shift of order S are the known spatial or angular 
Goos-Hanchen shifts for TM (p) or TE (s) polarization. In the '+/ O x ' block (and 
therefore also in the '+/ O x ' block) it is interesting to note that all the longitudinal 
shifts are equal to the general Goos-Hanchen shift d p or d s and that the transverse shift 
for all entries in this block are given by the same complex valued function u, the real 
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Table 2. Table giving the formalae for the symbols in Table 1. The formula is valid for 
both total internal and partial reflection, but the nature of the shift (angular, spatial 
or both) depends on the whether the formula yields a real, imaginary or generally 
complex number, and this is different for total internal and partial reflection. 



Symbol 




Formula 


Total internal reflection 


Partial reflection 


d p 




— / 

r 

p 

r 


imaginary 




real 


d s 




r 


imaginary 




real 


P 




l*>l 2 + |r s | 2 


imaginary 




real 


Q 


icot 


(r D +r s ) (r*+r* ) 
VO (|r p p + |r s |*) 


imaginary 




imaginary 


r 


cot 1 


n (r v -\-r s )(r*—r*) 
70 \r P \ 2 + \r s \ 2 


imaginary 




real 


u 




cot O r -4^ 


v* complex 




real 


V 




cot 


u* complex 




real 


a 




¥ 'p +¥ 's 

r p +r s 


Ima = — ip complex 




real 


b 




2 cos 6*o 


real 




real 


c 




r' -r' 

p s 

r p — r s 
^p+ir's 
r p +ir B 


Imc — —ip complex 




real 


e 




Im e — — ip complex 


9* 


complex 


f 


(1- 


hi) COt 0o Mr 


imaginary 


h* 


complex 


9 




t — ir 

' p 11 s 

r p —ir s 


Im g = — ip complex 


e* 


complex 


h 


(1- 


i) cote - p + r J 

! r p +ir s 


imaginary 


r 


complex 



and imaginary part of with contribute to the spatial and angular shifts of the circular 
and diagonal parts of this block respectively. More importantly for both circular and 
diagonal polarizations the transverse shifts are equal and opposite between the two 
orthogonal polarizations. 

The transverse part of the '+/0' block contains the spin-Hall effect of light [9], but 
from the symmetry of Table 1 follows that there is an analogous but reversed effects, for 
which the two incident circular polarizations experience shifts of equal magnitude but 
opposite sign when postselected by the same eigenpolarization. One noteworthy feature 
in the 'O x / O x ' block is the transverse shift for the ' x/x ' polarisation which is given 
by the longitudinal shift d s . In fact this is not solely a consequence of the algebraic 
approach, but depends on the form of the Fresnel reflection coefficients. 

Indeed more symmetries can be discovered if the functional form of the Fresnel 
coefficients and their derivatives are substituted for the algebraic expressions, but still 
without specifying if the reflection is total or partial. Taking this into account, there 
is a material independent shift of a component in the 'O/O' block, because b is a 
simple geometric factor b = 2cos6> - As this is real valued, the transverse shift in 
this block is purely spatial and its sign is given by the handedness of the incident 
polarization. Despite the fact that the other circular polarization has no transverse 
shift, the transverse shift for the total beam is not material independent, because the 
weightings in Eq. (2.5) depend on the Fresnel coefficients. 
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In total internal reflection f p and f s are unimodular (i.e. |r p | 2 = |r s | 2 = 1), and 
the ratios f' p /f p and f' s /f s are therefore purely imaginary. In this case d p and d s denote 
the purely spatial shift first discovered by Goos and Hanchen. Moreover, the function 
p = {r~*r' p + r*r' s )/2 in the row without analyser is also purely imaginary and there 
is no longitudinal angular shift for circular or diagonal polarisation either. Similarly, 
q = i(f p + f s )(f* + r*) cot 9 /2 and r = (f p + f s )(f* — ¥*) cot 9o/2 are purely imaginary 
and this means that there is no angular shift in the transverse direction either. This 
is a known result as the angular shift is related to the derivatives of the moduli of 
the reflection coefficients which in total internal reflection are zero [7] . In particular for 
incident circular polarisation the presence of a transverse spatial shift for the total beam 
is well known as the Imbert-Fedorov effect [3, 2]. 

In total internal reflection, u = v * = (f p + f s ) cot 0o/r p which is generally complex, 
and the scalar components can thus experience angular shifts while the total beam 
does not. Within the '+/ O x ' block incident circular polarisation thus gives rise to 
both transverse spatial and transverse angular shifts, but the angular shifts for TM and 
TE postselection average to zero because the total beam shows no transverse angular 
shift. Similarly, u and — v have equal but opposite real parts which account for equal 
and opposite transverse angular shifts between the TM and TE scalar components for 
diagonal incident polarization. In the 'O x / O x ' block, a,c,e,g all have the same 
imaginary part —ip, which means that all the longitudinal spatial shifts in this block 
are identical, which of course is reflected by the fact that the purely imaginary function 
p determines the longitudinal spatial shift for the total beam. 

In partial reflection from a dielectric interface the Fresnel coefficients are purely real 
and therefore d p and d s are as well, and represent only longitudinal angular shifts. Also, 
p is purely real in this case, and there is no longitudinal spatial shift of the total beam for 
any incident polarization. However, q is purely imaginary and the total beam therefore 
undergoes not only a longitudinal angular shift, but also a transverse spatial shift, which 
is of equal magnitude but opposite sign for incident left and right circular polarization. 
On the other hand r is purely real and the total beam for incident diagonal polarization 
experiences a general transverse and longitudinal angular shift but no spatial shift. In 
the '+/ O x ' block there is no longer a direct relation between u and v for partial 
reflection. Both functions are real and the TM and TE components of both circular and 
diagonal incident polarization experience a different transverse shift. In the 'O x / O x ' 
block a and c are real, while e = g* and generally complex. The presence of e and g 
in the 'O/x' block is the cause for the 'in plane spin separation of light' as discussed 
in [10]. Of course, because of the symmetry of the table a corresponding effect can 
be expected for a separation of diagonal polarisation components upon reflection of a 
circular polarised light beam and such an effect is evident in the presence of the same 
shifts in the 'x/O' block. 

A graphical representation of the shifts contained in this table is shown in Fig. 3 
for total internal reflection and in Fig. 4 for partial reflection. For both cases each row of 
Table 1 is compressed into two plots of overlapping circles, representing the spatial and 
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Figure 3. Tabic of the scalar component shifts in total internal reflection (n=2/3, 
Oo = 7r /4)- For each of the six incident polarization in the top row, the table shows the 
spatial and angular shifts for six corresponding polarization analyser settings, plotted 
on top of each other. The size of the rings is proportional to the reflected intensity 
of this component, normalised to the total reflected intensity. The shift is given in 
dimcnsionless unit and has to be multiplied by the appropriate measure (— 1/k and 
(<5 2 )) to give a physical quantity. The direction x and y are beam coordinates and refer 
to shifts with the plane (Goos-Hanchen) and orthogonal to it (Imbert-Fedorov) . The 
size of the circles is not indicative of the scaling between the shift and the spot size of 
the optical beam. The black dot and circle O indicates the position of the unshiftcd 
virtual reflected beam. The white dot and outline circle @ marks the position and 
intensity of the total reflected beam, without filtering out a scalar component using a 
polarisation analyser. 




Figure 4. Table of the scalar component shifts in external reflection (n=3/2, 
do = 7r/4). The table is constructed analogous to Fig. (3). Because in partial 
reflection the Fresnel coefficients re purely real, any imaginary units are introduced 
by the circular polarization which explains why spatial shifts only occur when either 
the polarizer or the analyser selects circular polarization. 
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angular component shifts for a given incident polarization, refractive index (n = 2/3 
for Fig. 3 and n = 3/2 for Fig. 4) and incidence angle of = ir/4 radians. The 
scalar component for each analyser polarisation follows the same colour coding as the 
incident polarisation and also as in Table 1. The position of the virtual reflected beam 
is indicated by a black circle and dot O, whereas the shift of the total beam is indicated 
by a white circle and dot @. The shifts are given as dimensionless quantities and have 
to be multiplied by — 1/k for the spatial shift and (5 2 ) for the angular shift. 

In addition to just quantifying the information given in Table 1 these figures also 
show the intensity of the reflected scalar encoded in the area of the each circle. The 
intensities are normalised with respect to the total reflected beam and should not be seen 
as indicative for the beam waist, which is usually so large that any spatial shifts are fully 
contained within its radius. Being able to read out the intensity of a scalar component 
helps to understanding how shifts of different magnitudes between an orthogonal pair 
can average to zero for the total shift. An example is the angular shift for circular 
polarisation in total internal reflection in Fig. 3, where the component with identical 
incident and analyser polarisations experiences a larger shifts than for cross-polarization, 
and yet the total shift is zero because of the different intensities in each component. This 
behaviour is fully in line with our understanding of beam shifts as realisations of weak 
values, where a lower intensity can go along with a larger value of the shift. 

5. Summary 

We have presented here a self-contained approach to optical beams shift phenomena 
which results in a simple linear algebra of 2 x 2 matrices, in which the connection between 
optical beams shifts and weak values is readily established. Moreover, on choosing three 
mutually incompatible bases for the incident and analyzer polarization ,we have been 
able to study the shifts of the centre of mass and the separate polarization components 
algebraically in a systemic manner, which resulted in Table 1. Using this table we have 
highlighted the existent of symmetries between input and output polarizations, which 
means that well known effects, such as the spin-Hall effect of light [9] and the 'in plane 
spin separation of light' [10] are accompanied by analogous but reverse effects. Further, 
we have found a material independent spatial, transverse shift, which can be observed 
for circular cross-polarization. As optical beams shifts can be used to gain information 
about the reflecting interface, the presence of a 'zero', which does not depend on the 
material is naturally very important. 

The connection of beam shifts and weak values is not new and has been pointed 
out before by Hosten and Kiwat [9] and phrased in a purely classical language by 
Aiello and Woerdman [8], but what hitherto has not been fully appreciated is, how 
exhaustive the interpretation of optical beam shifts phenomenon as weak values really 
is. The interpretation of the angular shift as imaginary part of the weak value in 
section 2 makes it obvious that the natural measure for the angular shift is the second 
moment of the spectrum. Using this as unit explains the difference in the angular 
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Goos-Hanchen (GH) shift between 3D vector beams [7] and 2D 'sheet' beams [26] and 
also accounts for the recently observed radial mode dependence for the angular GH 
and Imbert-Fedorov (IF) shifts for general Laguerre-Gaussian modes [27]. The vortex- 
induced shifts discussed in this paper can also be explained as weak values, but the full 
extent of the analogies between quantum weak measurements and optical beam shifts 
deserves a separate exposition elsewhere as it can discussed without reference to a given 
optical beam and material interface [15]. 
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